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We determine the evolving segregated or mixed morphology of charged-particle systems with long- 
range power-law interactions and overall charge neutrality that have been quenched to a low temper- 
ature. Segregated morphology systems are characterized by the size of uniformly charged domains, 
L(t), the particle separation within the domains, Zaa(£), the particle flux-density leaving the do- 
mains, J{t), the width of reaction zones between domains, W(t), the particle spacing within the 
reaction-zones, ^As(t), and the particle lifetime in the reaction-zones, r(t). Mixed morphology sys- 
tems are essentially one large reaction zone, with L ~ Iab ~ Iaa- By relating these quantities 
through the scaling behavior of particle fluxes and microscopic annihilation rates within reaction- 
zones, we determine the characteristic time-exponents of these quantities at late times. The mor- 
phology of the system, segregated or mixed, is also determined self-consistently. With this uni- 
fied approach, we consider systems with diffusion and/or long-range interactions, and with either 
uncorrelated or correlated high-temperature initial conditions. Finally, we discuss systems with 
particle-like topological defects and electronic systems in various substrate dimensions — including 
quantum-hall devices with skyrmions. 
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I. INTRODUCTION 

Reaction-diffusion systems without long-range interac- 
tions involve two diffusing species of particles that 
annihilate on contact, A+B — > 0. The temperature is as- 
sumed to be low enough that the reverse reaction can be 
ignored, so that the densities at late times are determined 
by the uncorrelated initial conditions. For equal initial 
mean-densities of the particles and anti-particles, there 
are two possible morphologies at late times (jj. The first, 
for spatial dimension d > 4, is a well mixed morphol- 
ogy in which mean-field dynamics applies to the particle 
density, dtp = —Cp 2 . This leads to p(t) ~ l/t at late 
times. The second, for d < 4, is a segregated morphology 
consisting of single-species domains of characteristic size 
Lit). The evolution is via diffusive currents feeding par- 
ticles from the domains into reaction zones where they 
are annihilated by anti-particles. The evolution of the 
mean domain density, p ~ i~ d / 4 , is easily demonstrated 
when the two species have equal diffusion constants, so 
that the density difference Ap = pa~Pb satisfies a linear 
diffusion equation Jl|,|| . It also applies more generally || . 

In the mixed morphology initial fluctuations decay 
faster than the mean-density and can be ignored. In the 
segregated morphology, long-wavelength initial fluctua- 
tions decay more slowly and asymptotically determine 
the domain structure. In that case, the profile of the 
domains and the structure of the reaction-zone can be 
understood through the particle fluxes and the annihila- 
tion kinetics of particle pairs [§-|7j , as well as on a more 
formal level j|. 

It has been recognized from the beginning [jjj that 
long-range interactions between oppositely-charged par- 
ticles and antiparticles change the evolution of a reaction- 



diffusion system — both by changing the initial charge- 
density fluctuations, and by changing the subsequent dy- 
namics. However, progress has been gradual due to the 
greater complexity of long-range interactions both in an- 
alytical models H and in computer simulations, and by 
the need to understand the force-free case first. Neverthe- 
less, progress in long-r ang e models has been made along 
a number of fronts |10j-|l8|| . For systems with long-range 
interactions but no diffusion, and with uncorrelated ini- 
tial conditions, Toyoki presented an analysis of a mixed 
morphology JTc| ] . Complementing that work, Ispolatov 
and Krapivsky considered segregated systems and 
simulated a variety of force-laws in d = 1. For both long- 
range interactions and diffusion, with uncorrelated initial 
conditions, a self-consistent model by Ginzburg et al. 
and a complementary scaling model 13 has captured the 
density evolution for n > d — 1, where n characterizes the 
force between particles via / ~ r~ n (see Eqn. (|l|) below). 
The n = d — 1 Coulombic case has also been treated by 
Ohtsuki [Q. High-temperature correlated initial condi- 
tions have also been discussed 

In this paper, we are interested in the late-time evo- 
lution of initially random distributions of charges in ho- 
mogeneous systems with overall charge neutrality. We 
use the scaling behavior, with respect to length-scale, of 
the initial charge fluctuations, the resulting large-scale 
currents, and the local annihilation rate of particles to 
determine the morphology as well as the time-exponents 
of the characteristic lengths. Our approach is based on 
the assumption that domain structures are characterized 
by only two lengths: their size L(t) and their characteris- 
tic particle separation, lAA(t), where the average density 
p ~ . This is comparable to the dynamical-scaling as- 
sumption of phase-ordering systems 19 and, with appro- 
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priate physical input, leads to a self-consistent descrip- 
tion of the system evolution. Between the domains, we 
allow for reaction zones of width W(t), within which both 
particle species are mixed with typical spacing Iab (t) and 
lifetime r. For segregated systems, charged currents are 
absorbed in the reaction zones. For mixed systems, the 
reaction-zone pervades the system. Our approach of bal- 
ancing currents with annihilation within reaction zones 
is in the spirit of the treatment of force-free systems by 
Redner and Leyvraz [Q; and for these systems our re- 
sults agree with the Renormalization-Group analysis of 
Lee and Cardy ||. 

We ignore correlations apart from the characteristic 
lengths Iaa, Iab, L, and W which determine the den- 
sities and sizes of the domains and reaction zones. As 
a result, our approach is insensitive to early-time dy- 
namics, unequal diffusion constants 0|, and motion of 
domain boundaries. All of these factors help determine 
the amplitude of the growth laws. Nonetheless, from the 
scaling properties of currents and lengths we can extract 
the asymptotic time-exponents. One benefit of our ap- 
proach is that its conclusions are robust to the details of 
the system, and that it provides a vividly physical pic- 
ture of the morphological evolution of charged systems 
with long-range interactions. 

We should emphasize what is new in this paper. We 
do not assume which morphology the system selects, but 
determine it self-consistently from the physics. We treat 
a large number of different cases with the same unified 
approach, including the well understood diffusion-only 
systems. We determine the reaction zone width and den- 
sity, and also discover non-trivial domain-edge profiles. 
For uncorrelated initial conditions, we discover several 
regimes where ballistic annihilation is the dominant an- 
nihilation mechanism. For high-temperature equilibrium 
initial conditions, we determine the appropriate initial 
charge-fluctuation spectrum and show how it modifies 
the subsequent evolution and morphology of the system. 

In the next section we introduce long-range inter- 
actions with overdamped dynamics between initially- 
random particles and anti-particles. For equilibrium 
initial conditions, we calculate the expected initial 
long-wavelength charge fluctuations, characterized by /i 
[Sec. II A, see also Fig. From the charge fluctuations, 
using a scaling form for the domain profiles near edges, 
we determine the scaling of charged currents [Sec. II B, 
see also Tab. Q and Fig. @|. We then consider various 



mechanisms of particle annihilation [Sec. II C , see also 
Figs. [5] and ||. We combine our results on currents and 
on annihilation, by balancing the currents with the an- 
nihilation rates in the reaction zones [Sec. II D|. This is 
sufficient to determine both the system morphology, and 
the time dependence of the domain size L and of the av- 
erage density, p [Fig. [| and Tab. We next consider 
the reaction zone in more detail, and determine its width 
W and particle spacing Iab 1 &s well as the typical lifetime 
r within the reaction zone [Sec. HE, see also Tab. III|. 
Finally we discuss our results, including implications for 



coarse-grained treatments [Sec. [IF]. 

We discuss the previous reaction-diffusion literature 
[Sec. Ill], then discuss applying our results to electri- 
cally charged (n — 2) systems in various substrate di- 
mensions [Sec. IV]. We emphasize that for d < 3 novel 
decay laws and segregated morphologies are found. We 
then discuss the application of our results to quenched 
phase-ordering systems with point-like topological de- 
fects [Sec.0. We can also extend our approach to include 
Levy super-diffusive systems, as well as to sub-diffusive 
systems [Sec. |Vlfl . 

The effect of short-range cutoffs on the power-law in- 
teractions can easily be treated [Sec. IVIJ]. We numer- 



ically explore our results on domain profiles [Sec. VIII 
see also Fig. |(|. Finally, we conclude [Sec. IX] 



Throughout this paper we concentrate on exponents, 
or scaling dependencies, of various quantities in the late- 
time limit. Inequalities apply to exponents, so that a 
process is dominant if it is asymptotically largest as t — > 
00. We denote the mechanisms with a subscript D, F, or 
B for diffusive, local force driven, or long-range ballistic 
processes, respectively. 



II. LONG-RANGE INTERACTIONS 

We work in the overdamped limit, in which particle ve- 
locity equals the applied force times a constant mobility, 
n. The pairwise forces are 



fij = {n-l)Cq iqj /r%, 



(1) 



between particles with charges {q{\ and pairwise separa- 
tions {rij}. This corresponds to a pairwise interaction 



energy E i: 



Cq,q 3 



with a logarithm at n = 1. 



Thus the particle velocity is 



, fijUj + 



(2) 



where we have added a random uncorrelated noise for 
diffusive motion with diffusion constant D, where (4>i(t) ■ 
4>j(t')) — DSijS(t — t'), When oppositely charged parti- 
cles get within a fixed capture radius r c , they annihilate 
instantaneously pit]. We are interested in the behavior 
of the system at late times when distinct length-scales 
are well-separated. 



A. Initial Conditions 

Most studies of reaction diffusion with long-range 
interactions have focused on the case of random un- 
correlated initial conditions, where the particles are 
randomly placed with a local Poisson distribution 
[|l^-0,|6Ul|Jl^,|l|] . Experimentally, it is more natural 
to quench a system of charged particles from a high- 
temperature state in which they are in thermal equilib- 
rium [fpJil. 
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Charge fluctuations may be usefully characterized by 
the typical excess charge density at scale L, 

Sp-L-f*. (3) 

This can be thought of as the excess charge density in a 
region of size L after coarse-graining to that scale. It is 
also related to the magnitude of the Fourier-component 
Pi /l of the charge-density by 

\ Pl/L \~L d / 2 -», (4) 

as obtained by summing up 0(i^,/ 'L d ) uncorrelated con- 
tributions of size L d Sp(L) and normalizing out the con- 
stant contribution due to the system size, L^. As a result 
we can obtain p directly from pk- 

We are only interested in the scale dependence of initial 
charge fluctuations that survive for long times — i.e. for 
fluctuations at large scales and correspondingly small fc. 
Hence we consider a continuum charge-density p(r) with 
energy 

E = C/2 J d d rd d r'p{r)p{r')/\r-r'\ n - 1 

From efc = J d d xe lk ' x x 1 ~ n , we obtain = 
7T d / 2 (fc/2) n - 1 - d r[(d+l-n)/2]/r[(n-l)/2] for 1 < 
n < d + 1. For n > d + 1 we must introduce a UV 
cutoff, the inverse particle separation p2]| , which deter- 
mines a leading ej, = const small-fc behavior. For n = 1 
the interaction in Eqn. (|^) should be logarithmic, and we 
get tk = 7r d / 2 2 d ~ 1 k~ d T(d/2). Imposing equipartition on 
Eqn. (||) with temperature T, we get (pkP-k) ~ ksT/ek- 
This derivation reproduces the small k behavior of a 
more complete variational approach, see, e.g. |23f| . En- 
tropy factors do not affect the leading long-wavelength 
fluctuations, so that the integrand in Eqn. (H) is correct 
for small k. This means that pk is Gaussian distributed 

for k —>■ 0, and (\pk\) ~ ^ 2 '■ Comparing with Eqn. (^) 
we obtain 

/ (2d + l-n)/2 I<n<<2+1 , R , 

M_ \d/2 n>d+l. [b) 

This is illustrated in Fig. [j]. Note that d/2 < p < 
d, so that long-range interactions always suppress ini- 
tial charge-fluctuations. The maximal suppression is 
achieved at n = 1 when p = d. [For n < 1 we expect 
higher-point correlations to be significant.] For sharp 
enough interactions, with n > d+1, we recover uncorre- 
lated initial conditions with p — d/2. 

The charge excess at scale L, as discussed above, is 
quite different from the charge excess within a sharp 
boundary of scale L, for example inside a sphere of radius 
L. The later has been proposed as a measure of charge 
fluctuations, particularly in Coulombic systems in gen- 
eral dimensions, with n = d — 1, where a Gauss's law 



applies (see, e. g., p | JT5|Jl^ |) and in systems with topologi- 
cal defects [|24|-|26| where similar integral identities apply 
to the topological charge density. At high temperatures 
the charge inside a given closed surface is proportional to 
the square-root of the surface obtained from in- 

tegrating the appropriate random high-temperature field 
over the surface pq| . Taken literally this would imply 
that p — (d + l)/2 |l5|| , i.e. larger charge fluctuations 
than indicated in Eqn. (^|) for Coulombic interactions. 
However a sharp surface picks up charge fluctuations at 
short scales in addition to the desired large-scale fluctu- 
ations. Indeed, the Coulombic or high-temperature na- 
ture of the system is moot — arbitrary sharp surfaces 
within a system with microscopic charge heterogeneities 
(e.g. atoms or thermal fluctuations) pick up charge fluc- 
tuations proportional to the square-root of the surface 
area. The useful and appropriate charge fluctuations at 
scale L are given by Eqn. (^), which are not mixed with 
charge-fluctuations at short-scales [^7j . 

B. Domain Profile and Currents 

Charge fluctuations, coarse-grained to scales much 
larger than typical particle separations, decay via 
charged currents. These currents can be diffusive, or can 
be driven by the long-range interactions. We first con- 
sider segregated systems, where the domain size, L(t), 
sets the scale of surviving charge fluctuations. The av- 
erage density within a domain, p ~ is simply pro- 
portional to the initial fluctuations at the domain scale, 
5p(L), so that m 

L ~ l d H (7) 

Charge fluctuations at scales much larger than L cannot 
have relaxed since charge transport at larger scales is cut- 
off by the domain structure. We now assume p8[ that the 
domain profile is only determined by the average density, 
p, and size, L, of domains. For example, domains have a 
typical density profile 

p(r)=l AA f(r/L), 

~ll d A (r/L) a , (8) 

where r is measured from the edge of the domain. The 
second equation holds near a domain edge, with r <C L. 
We also require r ^> W so that the profile is probed well 
away from the reaction zone. This scaling form has been 
proposed by Leyvraz and Redner in diffusive systems |Q , 
and was numerically confirmed in d = 1. 

The assumption of an invariant scaled domain mor- 
phology is powerful, and is sufficient to determine the 
coarse-grained current density J near the domain edge. 
The flux must have a dominant non-zero constant contri- 
bution near the domain edge for r -C L arising from the 
evolution of domain density, since no annihilation takes 
place in the domain interior. Consider the net charge 
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of a domain, Q ~ L d p ~ L d ~^. It implies a non-zero 
net flux density J ~ Q/ L^ 1 ~ L 1 ^^ jt near the domain 
edge. The exponent characterizing the domain profile, a, 
is constrained to allow the dominant J to be finite but 
non-zero for r/L — ► + . 

If the dominant current is diffusive, then Jd ~ so 
that 



(9) 



Imposing a constant diffusive current condition at the 
domain edge implies a linear profile, with a — 1. This 
agrees with studies of domain profiles in diffusive systems 



If the dominant current is driven by long-range forces 
then it is given by the coarse-grained field times the local 
charge density, Jf ~ p{r)F(r). The field, F{r), at a 
distance r away from the domain edge due to the charge 
distribution given by Eqn. (|J) is 



F(r) 



Iaa(t) 



d d x[p(r +x)- p(r - x)]/x r ' 



(10) 



We have used the domain scale as the long-distance 
cutoff, and the local interparticle spacing, Iaa{t) ~ 
p(f)~ l / d ', as the short-distance cutoff. We also restrict 
the angular integral to be close to the normal direction 
from the interface, which retains the scaling behavior of 
F(r) without requiring detailed information about the 
domain shape. 

For d > n, charges at distances of order L determines 
F, so that F ~ pL d - n ~ L<l-n-^_ For a constant flux 
Jf at small r, we must have p(r) ~ const, (i.e. a = 0) 
so that J F ~ p 2 L d ~ n ~ £<*-n-2/i_ when d > 2n, with 
uncorrelated initial conditions (p — d/2), F increases 
with the upper-cutoff of the integral in Eqn. ([To]) , so that 
we should use the system-size, L^, rather than L. For 
this case, the system size enters the dynamics, and the 
thermodynamic limit does not exist juj[[l) . 

For d < n, the local charge distribution dominates the 
F(r) integral. We can expand the density for x <C r and 



find F(r) ~ p(r 



,(n-l)/d, 



for n > d + 1 . For d < n < 



d + 1, on the other hand, the integral is dominated by 
scales around r, and we find F(r) ~ r a+d ~ n /[l\ A L a ]. 
Insisting on Jf approaches a non-zero constant near the 
domain edge, we obtain Table [jj These results apply 
far from the reaction zone, but otherwise close to the 
domain edge: W <r <C L. [We determine the reaction 
zone width in Sec. HE.] We see that currents dominantly 
driven by power-law interactions lead to non-trivial do- 
main profiles, with § < a < \/2, in dramatic contrast 
to the diffusive case where a = 1. It is also interest- 
ing that the current, Jp, has the long-range form for 
d < n < d + 1, even though a ^ 0. 

When both long-range forces and diffusive effects are 
present, then we must compare Jf and Jd, and identify 
which is larger at late times. We summarize the results 
in Fig. |^ for random uncorrelated initial conditions. For 



equilibrated high-temperature initial conditions Jd is al- 
ways asymptotically larger at late times. 



C. Particle Annihilation 

In a well mixed region of the system, where the typi- 
cal spacing between oppositely charged particles is Iab, 
what is the scaling of the particle lifetime t(Iab)^ There 
are three annihilation mechanisms: diffusive annihila- 
tion (td), local interaction-driven annihilation (rp), and 
ballistic annihilation (tb). We determine their scaling 
dependence on Iab and hence identify the dominant 
mechanism at late times. Our essentially microscopic 
approach also provides insight int o the applicability of 
coarse-grained treatments (see Sec. [IF). 

In the force-free case, particles move diffusively with 
diffusion constant D, and annihilate with oppositely 
charged particles when they approach within a fixed dis- 
tance r c . In d < 2, trajectories are space filling and 
td ~ Iab/D — the time it takes for a particle to diffuse 
Iab- For d > 2, td ~ Iab I '(Dr d ~ 2 ) , since each parti- 
cle must explore the characteristic volume per particle to 
find an anti-particle to annihilate. Wc have 



td 



I 2 

AB 
id 
AB 



d<2, 
d>2. 



(11) 



In the diffusion-free case, considering only local inter- 
actions, / ~ r~ n , between two particles initially sepa- 
rated by Iab, the annihilation time 



TF 



jn+1 
AB ■ 



(12) 



For many particles in a region, the same result follows 
from the scaling of the velocities in Eqn. (^). 

With both diffusion and local interactions, diffusion 
(td) dominates the annihilation time for n > 1, while for 
n < 1 the force (rp) does. This follows directly from the 
particle dynamics, Eqn. (^). Rescale all distances by Iab, 
and rescale time by l\ B , so that diffusion is unchanged in 
the scaled coordinates as Iab increases. Scaled velocities 
due to the force are then multiplied by l\' B n . As a result, 
forces do not asymptotically contribute for n > 1, while 
forces dominate for n < 1. For d > 2, this leads to the 
initially counter-intuitive result that diffusion dominates 
for 1 < n < d — 1, even though td ^> Tf. This is a 
well-known result for n — 2 (g9|. Essentially, the com- 
petition between diffusion and the attractive interaction 
is along the separation vector between two particles, and 
hence is always one-dimensional in character. Indeed, in 
d = 1 the faster annihilation mechanism dominates, and 
the marginal value is n = 1. We have confirmed these 
predictions for various n in d < 3 by placing a particle 
and an anti-particle in a periodic box, and plotting the 
average annihilation time as a function of the box size. 
In Fig. ||, we show our results for d = 3. For n > 1, 
the long-range force merely changes the effective capture 
size, leaving td ~ Iab- 
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When an applied or non-local force F is present, we 
must also consider ballistic annihilation. With a velocity 
proportional to F, and in a time Tb, particles sweep out a 
volume proportional to tbFI% , where Z* is the radius of 
the effective capture cross-section. Equating this to the 
typical volume l d AB per particle, we obtain the typical 
ballistic annihilation time 



tb 



id ,1- 

L AB L * 



l /F. 



(13) 



When no noise is present, a particle will be captured by 
local interactions only at separations less than Z* where 
local interactions are as large as F, so that F ~ / ~ 1/7" 
or can use the typical forces F(r) from 

Table [|, but must evaluat e the m at the reaction-zone 
width W (see below, Sec. HE) — since the reaction- 
zone itself is neutral at a coarse-grained scale. Because 
F depends on the system morphology, this requires a self- 
consistent solution [30t| . The results are simple: for n > d 
we find Z* > Iab and the dominant time-scale is Tp , while 
for n < d we find I* < I^b and so the dominant time-scale 
is tb- This corresponds to the naive phase-space result 
from Eqn. (|l^) that charges from far away dominate local 
forces, and hence tb dominates, only when n < d. 

When diffusive noise is present, we first ignore the lo- 
cal interactions and only consider an applied field F. For 
d > 2, both the random walk and the ballistic force sweep 
out volume proportional to time. Because F decreases 
with time, td always asymptotically dominates. How- 
ever, for d < 2 random walks recur, and the ballistic drift 
can enhance the volume covered by the random walk by 
suppressing the recurrence. The rate of volume swept 
out by the drifting particle, is a sphere of radius Z* every 
At ~ U/F ~ ll/D. This implies Z* ~ l/F, and leads to 

id J?d—2 
TB ~ l AB f 

When both noise and local interactions are present, 
the local capture cross-section Z* is determined by the 
dominant non-ballistic process — i.e. Td for n > 1 and 
tf for n < 1 . The shortest annihilation time for uncorre- 
cted initial conditions is given in Fig. ^|. For equilibrated 
initial conditions, td always dominates for n > 1. 

Our treatment of particle annihilation is essentially mi- 
croscopic rather than coarse-grained, since we have built 
the particle separation, Iab, directly into the annihilation 
times. We have derived our results by considering parti- 
cle pairs, though we have included some multi-particle ef- 
fects by never allowing particles to "escape" further than 
Iab from an anti-particle. We apply the results in mixed 
regions of the system such as reaction zones. 



D. Domain Morphology: Segregated or Mixed 

We assume the system is described by one of two 
morphologies pi] , depending on whether coarse-grained 
charge fluctuations are comparable to or much less than 
the mean particle density at late times. The former case 



describes a segregated morphology with domains of parti- 
cles separate from domains of antiparticles. In the second 
case, there is a mixed morphology with no clearly defined 
domains. We can use the dominant flux J and fastest an- 
nihilation time t to see which morphology is consistent. 
The system turns out to have a unique consistent mor- 
phology: either mixed or segregated. 

First consider a segregated morphology. The system 
has domains of scale L separated by reaction zones of 
width W. Within the zones there is a typical particle 
spacing Iab- We ignore correlations or structure within 
the zones [p^ . Our self-consistency constraints are that 
Iab ^ Ua, and that W < L. We impose the former 
because annihilation takes place in the reaction zone but 
not in the domain bulk, so the density in the reaction 
zone should be smaller as a result. We impose the lat- 
ter since if it were not the case, with W ^> L, the 
system would be effectively all reaction zone and of a 
mixed morphology. These constrain the maximum rate 
that reaction zones can "process" incoming particles: the 
average density of particles that are in reaction zones, 
WL d ~ x /(l d AB L d ), is at most l/l AA . Conversely, the max- 
imum annihilation rate is 1/t(I aa ), since r(Z) is mono- 
tonic. Combining these, the maximum annihilation rate 
from reaction zones is dtpmax ~ ^1/[I aa t(Iaa)]- For 
segregated systems, the actual rate of particles annihi- 
lating per unit volume is determined by the currents en- 
tering the reaction zones: dtpseg ~ — JL d ~ x /L d ~ —J/L. 

If the maximum reaction rate dtpmax asymptotically 
dominates the actual rate dtp se g, then a domain structure 
is consistent. Indeed, since dtpmax is the rate of density 
decrease in a mixed morphology, the inequality indicates 
that the background density decay quickly with respect 
to the charge fluctuations — i.e. that the segregated 
structure occurs when it is consistent. If the maximum 
reaction rate is less than the particle flux, then segregated 
domains cannot be sustained, and the mixed morphology 
should result. 

For the segregated morphology, we use dtp ~ — J/L ~ 
—p/t and L ~ l AA to extract the domain scale and par- 
ticle density within the domain. 

For the mixed morphology, the flux does not drive an- 
nihilation of the mean density. Rather every particle is 
effectively in a reaction zone and has a characteristic life- 
time r. Comparing this to the scaling of the density evo- 
lution, dtp ~ —p/t ~ ~P/t, indicates that t(Iaa) ~ t. 

The resulting growth law regimes are shown in Fig. |^ 
for systems with either uncorrelated or high-temperature 
equilibrated initial conditions, and with either force-only 
or force-and- noise dynamics. The exponents are sum- 
marized in Tab. H. The noise-only case reproduces the 
Ovchinnikov-Zeldovich-Toussaint-Wilczek result M , and 
is given in Fig. || b) by the short-range (n — > oo) limit 
of the force-and-noise dynamics with uncorrelated initial 
conditions. 

We can now consider the consistency of the mixed mor- 
phology. The system has no domain structure, and has 
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local charge separation scales Iab ~ Iaa- Coarse-grained 
charge fluctuations, above a scale X(t), remain from the 
initial conditions. The mixed morphology is consistent if 
the charge density of the remaining charge fluctuations 
Sp(X) <~ A _Al is much less than the mean particle density 

,-d 

P ~ AA' 

Given the scaling of net currents J(X), the charge fluc- 
tuations evolve with d t Sp(X) ~ Sp/t ~ JX d - 1 /X d ~ 
J IX. For diffusive currents, J D {X) ~ Sp(X)/X ~ 
and we find that X ~ i 1 / 2 . For force driven 
currents, F(X) ~ x d - n -» from Eqn. @ using X as 
the upper and lower cutoff. The force drives the net 
charge density Sp(X), and yields a net current Jp{X) ~ 
leading to X - ^/(M+n+l-d), Using tne re _ 

suits for Iaa from above, we confirm that Sp(X) -C p 
within the mixed regions of Fig. |^. 

Our exponents apply at regime boundaries, where they 
are continuous. However our approach does not deter- 
mine any logarithmic factors. [Indeed, logarithmic fac- 
tors are expected at n = d, whenever Jp or tb domi- 
nates, due to the logarithmic divergence of the integral 
for F, Eqn. ([h]).] At a boundary between mixed and 
segregated regimes, any logarithms present determine 
the dominant morphology through the self-consistent ap- 
proach described above. Without logarithms, the mixed 
morphology applies on the boundary, since W ~ L there. 
The width W characterizes the mean-distance to annihi- 
lation for particles, and thus there is a finite density of 
anti-particles any finite multiple of W into a domain. 
When W ~ L, there is a finite density of anti-particles 
arbitrarily deep within a domain — i.e. the system is 
mixed. This is the case in the diffusion only case at 
d = 4, which mixes M. 



E. The Reaction Zone 

For reaction-diffusion systems with a segregated mor- 
phology, much progress has been made on the structure 
and evolution of the reaction zones between domains 
§-§. With Ion g-range interactions, we develop a similar 
approach that balances the dominant flux into the re- 
action zone with the dominant annihilation mechanism 
within the reaction zone. This balance determines the 
reaction zone width, W(t), and the interparticle spacing 
within the reaction zone, Iab (t) ■ To simplify our discus- 
sion, we take the density of both species to be uniform 
throughout the reaction zone p2| . 

When a particle enters the reaction zone, it travels a 
typical distance W before it annihilates — the "annihila- 
tion mean-free path" . For a segregated system, W must 
characterize the width of the reaction zone. A reaction- 
zone that grows less rapidly than W allows almost all 
particles to pass through unannihilatcd. On the other 
hand, a reaction-zone that grows more rapidly than W 
would provide infinite number of annihilation mean free- 
paths as t — > oo, and would not have a mix of particles 



and anti-particles at the far edge of the reaction-zone. 
Since slower or faster growth is not self-consistent, W{t) 
characterizes the typical reaction zone width. 

In the annihilation time t, a particle diffuses a typ- 
ical distance Wd ~ (Dt) 1 ^ 2 , or it ballistically moves 
Wb ~ Ft under an applied field F, where F is given 
by Table |. In any case, a particle must move at least 
the interparticle spacing to annihilate, Wf ~ Iab- The 
largest of these widths describes how far a particle trav- 
els before annihilation, and so characterizes the width of 
the reaction-zone. 

The reaction lifetime, r, the particle spacing in the 
reaction-zone, Iab, and the zone width, W, are de- 
termined self-consistently. Given the dominant flux J, 
we equate the overall flux density into reaction zones 
Nfiux ~ JU 1 ^ 1 /L d with the rate of annihilation within 
the reaction zones WL d ~ l /[L d l d AB T(l AB )]. We choose 
the largest W for the given r, and the fastest r for the 
given Iab- These dominate the width and annihilation, 
respectively. This provides a self-consistent solution for 
r, Iab, and W. 

For diffusion-only systems in d < 2, our argument is 
essentially that of Leyvraz and Redner f3]. We recover 
their results of W ~ Iab ~ t 3 ^ 8 , an d t 1 / 3 in d = 1 and 2, 
respectively. For d > 2, however, we have W ~ y / r5" ~ 

l AB since random-walks are no longer space-filling, and 



d = 3, and 



we obtain 1 A b ~ i 5/18 and W ~ t 5/1 
Iab - Ua ~ 1/4 and W ~ L ~ t 1 / 2 in d = 4. Our 
results coincide with the RG scaling results of Lee and 
Cardy ||. We also recover their scaling relations with 
respect to currents, whenever td dominates annihilation. 
In particular, with dtp ~ J/L ~ W/(LI ab td), t b from 
Eqn. and Wd ~ \[td, we have Wd ~ J -1 / 3 and 

Iab - J- 2 '^ for d > 2 and W D - Iab - J-^^+V for 
d < 2. These relations hold independently of the process 
that dominates the currents. 

Even within the mixed morphology, it is interesting to 
consider the distance W a particle travels in a lifetime 
t. This distance W(t) is the same scale as the remain- 
ing charge fluctuations, X(t), as discussed before in the 
previous section. This is reasonable, since charge fluc- 
tuations at scales much larger than W cannot be flat- 
tened out by charge motion. We summarize our results 
in Tab. fu. 



F. Discussion 

In all cases the density decays more slowly as n in- 
creases and the potential becomes sharper and hence 
shorter ranged. If noise is present, then diffusive pro- 
cesses eventually dominate as n increases. As the spatial 
dimension d is increased for any n, then a mixed mor- 
phology is eventually reached. For generalized Coulombic 
systems, with n = d— 1, the density decays as p ~ l/t in 
all combinations of high-temperature or uncorrelated ini- 
tial conditions and force and/or diffusive processes. They 



G 



are always of mixed morphology. 

With equilibrated initial conditions, diffusive processes 
dominate independently of n when they are present. This 
is reasonable, since the high-temperature initial condi- 
tions balances the interactions with temperature. Tem- 
perature then becomes more relevant as it is reduced dur- 
ing the quench. 

For systems without diffusive processes, the growth 
regimes and processes for uncorrelated and high- 
temperature initial conditions, in Fig. [b| a) and d) respec- 
tively, correspond. The various growth exponents differ 
only due to the different initial charge fluctuations, char- 
acterized by p. On the other hand, with both diffusive 
and long-range processes the regimes do not correspond, 
between Fig. [| b) and e), because the competition be- 
tween force and diffusive processes depends on the charge 
fluctuations through p. 

Exponents are continuous at regime boundaries for 
Iaa, Iab, L, r, and W. For segregated systems L > 
W > Iab <L Iaa, where the last inequality holds since 
reactions should decrease the particle density. We also 
check that r < t in all cases. In d = 1, we check W ~ Iab, 
as expected since the reaction zone is precisely one AB 
pair. At the border between mixed and segregated mor- 
phologies, the reaction zone is maximal, i.e. W ~ L and 
Iab ~ Iaa- 

The particle density at the edge of a domain, at r ~ W 
in Eqn. (|J), scales the same as the particle density within 
the reaction zone, l A d B . i.e. that l~ A d B ~ l A d A {W/L) a . In- 
teresting special cases occur where Jf dominates the flux 
but Wd dominates the reaction zone width (regions [ii 4 ] 
and [ii 5 ])- In these special cases, the singular domain 
profile, with a < I , leads to a diverging diffusion flux as 
the domain edge is approached, which must eventually 
dominate at some distance Lx from the domain edge. In 
particular, Lx is determined by p{Lx)/Lx ~ Jf- Using 
the domain profile, Eqn. (||), we have L~ >1 ~ ol L x ~ ~ J F - 
If Lx <C W, then the cross-over is pre-empted by the 
reaction zone and we do not expect to observe it [^3| . 
However, precisely when Jp and Wd dominate we find 
that Wd -C Lx < L, so that a linear transition regime 
is expected for Wd <C r <C Lx- The intermediate lin- 
ear regime, when it occurs, does not change the density 
evolution or reaction-zone structure, since Jf still char- 
acterizes the flux. 

In reaction-diffusion systems, d c is the critical dimen- 
sion above which a coarse-grained reaction rate XpaPb 
applies, and d u can be defined as the dimension above 
which we can neglect inhomogeneities. Below d c the re- 
action term is not given by p 2 in a coarse-grained de- 
scription, and local density fluctuations must be taken 
into account |^,^|. Above d u we can ignore spatial gra- 
dients and have dtp ~ — p 2 and hence p ~ 1/t (sec jjj). 
Both of these effects stem from the diffusive annihilation 
mechanism. The local annihilation rate p/td is propor- 
tional to p 2 only for d > 2 where td ~ p^ 1 ■ This sets 
d c = 2. The mixed state is found for d > 4, so that 



d u = 4. 

When only long-range interactions are included, these 
definitions are not as useful since neither tf nor tb are 
in general proportional to p^ 1 . As an example, when 
only force-driven evolution is included, as in Fig. || a) 
and d), the density in the mixed state has a d-dependent 
exponent, see [Hi] and [vii]. Even the long- wavelength 
charge fluctuations remain relevant to drive the dominant 
ballistic annihilation. We could describe this as d c = 
d u = co. 

When diffusive and long-range processes are included, 
however, we see from Tab. Ill that the diffusive process 
t d ~ Iab dominates above d c = 2, so that the critical 
dimension is unchanged from the diffusion-only case. In 
mixed morphologies diffusive annihilation always domi- 
nates, so that d u equals the dimension where the mixed 
morphology starts. For uncorrelated initial conditions we 
find d u — n + 1 for n < 3, while d u = 4 above. For high- 
temperature initial conditions d u = (n + 3)/2 for n < 5, 
while d u — 4 above. Since diffusive processes dominate 
all regimes in Fig. ^| e) , it is only through the suppression 
of initial charge fluctuations that the long-range interac- 
tions modify d u - 

We have assumed that the domains of the segregated 
morphology are characterized only by a size L and char- 
acteristic particle separation Iaa, and that the reaction 
zones are similarly characterized by a width W and a 
particle separation Iab- This leads to the exponents 
summarized in the tables and figures. It remains pos- 
sible that other lengths enter into the asymptotic evo- 
lution, for instance through fluctuations in the shape of 
domain boundaries (see, e.g., ||). Our results are most 
robust for n < d, where only the domain scale L en- 
ters in the calculation of the forces and currents. Be- 
cause of that, we obtain the same result for the evolving 
density in a segregated morphology after a coarse grain- 
ing to O(L), and any new but shorter lengths would not 
change this result. It remains possible, however, that 
additional lengths change the structure of the reaction 
zone, and hence change the boundary between segregated 
and mixed morphologies. Similarly, higher point correla- 
tion effects may lead to more intermediate morphologies 
within what we label the mixed regime. For n > d addi- 
tional lengths would affect the evolution of the density, 
since short scales enter into the calculation of the char- 
acteristic forces and currents. 

Spatial fluctuations of coarse-grained quantities have 
also been neglected in our treatment. We assume that 
fluctuation effects are not strong enough to change our 
results for exponents. Another way to view fluctuations 
is in terms of the distribution of various quantities. We 
effectively assume that distributions have negligible tails. 

We have also not included the motion of domain 
boundaries, either from local heterogeneities in domain 
density, or from differing mobilities of the particle species. 
Domain motion should not affect the scaling exponents 
if the boundaries move slowly enough for the domain 
profile a to maintain itself. A simple check is reas- 
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suring. The flux required to move domain interfaces, 
J motion ~ pL ~ pL/t, has the same scaling behavior as 
the characteristic flux, J ~ Ldtp. Hence domain wall 
motion should introduce no new scales. 



III. PREVIOUS LONG-RANGE WORK 

Some of the earliest work on charged A + B — > sys- 
tems with long-range interactions was by Toyoki |Io[ . He 
treated uncorrelated initial conditions and force-only evo- 
lution — corresponding to Fig. || a). He presented the 
mean field analysis of the mixed morphology, our region 
[Hi]. By considering the mean-square force on a particle, 
he recovered the system-size dependence for n < d/2. He 
also initiated numerical studies of d = 2 systems. How- 
ever he used a short-scale cutoff pro port ional to the aver- 
age particle separation Iaa (see Sec. VII), which hampers 
interpretation of his results. Jang et al. | fl6| numerically 
studied d = 2 systems with n = 1 and obtained p ~ £-°- 55 
using a large noise amplitude, and p ~ £-0.90 usm g a 
large force amplitude. This is consistent with our results 
of p ~ i -1 / 2 and p ~ for the pure noise or pure force 
systems, respectively. 

Ispolatov and Krapivsky fllj focused on force-only evo- 
lution in d = 1, with a segregated morphology. They ob- 
tained results consistent with ours for n < d, where their 
assumption of a constant domain profile, with a = 0, is 
valid (see Table | and Sec. |yillp . 

Ginzburg, Radzihovsky, and Clark j^] treat coarse- 
grained hydrodynamic systems with uncorrelated initial 
conditions and with both diffusion and long-range forces, 
for n > d — 1. Our particle-based treatment can be seen 
as complementary to their work. We obtain the same 
results for density evolution, but we also obtain details 
about the reaction zone — including W , Iab, and r. We 
also identify the mixed or segregated morphology of the 
system. Furthermore, the system-size dependence in the 
force for n < d/2 fl(],[ll| is missed in their treatment. 
It would be interesting to extend their approach to in- 
clude the particle nature of the charges, possibly with 
new phcnomcnological annihilation terms. Burlatsky et 
al. |13] present a simple scaling model that obtains the 
same results as the self-consistent approach of Ginzburg 
et al. lil. 



IV. ELECTRONIC SYSTEMS 

Asymptotic non-geminate pair-recombination in clean 
crystalline semiconductor systems should be described 
by our approach with n = 2 and equilibrated high- 
temperature initial conditions. In d = 3, we find p ~ t^ 1 
with a mixed morphology. This result applies for every 
combination of uncorrelated or high-temperature initial 
conditions and force and/or diffusion. This is because all 
annihilation mechanisms have the same scaling, r ~ l AB . 



In contrast, for d < 3 the annihilation mechanisms dif- 
fer and the morphologies are all segregated. The structure 
of the evolving charge-fluctuations is thus much richer. 
Specifically, in d = 1, and 2 we expect regime [um'] to 
apply asymptotically with p ~ i" 1 / 4 and t~ 3 / 4 respec- 
tively. The decay rate is dramatically slowed due to the 
segregated morphology. 

Clean d = 2 systems exist in quantum hall effect 
(QHE) devices. QHE devices with skyrmion charge exci- 
tations (see, e.g., Q) may be particularly good systems 
to study these effects, due to their low mobilities and slow 
dynamics combined with sensitive time-resolved probes 
of their particle density. We will develop this in more 
detail in a separate publication (35), P avm g particular 
attention to the scaling of the amplitudes and the result- 
ing pre-asymptotic cross-overs in the evolution. 



V. PARTICLE-LIKE TOPOLOGICAL DEFECTS 

Point-like topological defects, e.g. hedgehogs in d = 3 
or vortices in d — 2, are found in liquid crystals and 
in vector O(N) systems in d = N dimensions. Consid- 
ered pairwise, these defects have power-law interactions. 
Indeed, early theoretical work [[l0| on these systems was 
based on the interactions of point-like topological defects. 
However the evolution of systems with point-like topolog- 
ical defects is different than the dynamics of Sec. || (see 
however p6[). The order parameter field that supports 
topological defects provides a scale-dependent mobility to 
particle motion. Additionally, the mobility depends on 
the local environment, i.e. the interaction between de- 
fects is not a two-point particle-particle interaction. We 
may, however, take that as a first approximation. For 
d = 2 XY or nematic systems, the mobility scales loga- 
rithmically with the particle velocity, 77 (Z ) ~ l/\n(dl/dt) 
p9| . For d = 3 0(3) or nematic systems, the mobility 
of hedgehogs scales as the inverse separation, r)(l) ~ l/l 
p7f . In general, for O(N) models in N > 3 dimensions, 
j](l) l 2 ~ N High-temperature equilibrated initial 

conditions are appropriate for quenches in physical sys- 
tems. For point defects in d > 3, the effective interaction 
is linear [[[(J (n — 0), which has not been treated in this 
paper. However, d = 2 systems are within our purview. 

For the d = 2 XY model, interactions between defects 
are logarithmic (n — 1). Ignoring the logarithmic mo- 
bility, and including equilibrated initial conditions p7| , 
we have a mixed morphology with all lengths scaling 
the same, Ua ~ Iab - L - W ~ t 1 ' 2 . This is found 
whether or not diffusive processes are included, and in- 
deed To ~ tb ~ Tp ~ l\ B and Jjj ~ Jp ~ L~ 3 for 
this system. The single length scale and the similarities 
between no noise quenches and quenches with diffusion 
matches the phase ordering results of the d = 2 XY model 
©|| . 

Two phase-ordering systems that do exhibit strong 
scaling violations are the d — I XY model |39| and the 
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d = 2 0(3) model [Q. Both of these systems support 
non-singular topological textures that have particle-like 
aspects. However topological textures have an intrinsic 
length-scale that evolves in time. Significant extensions 
to our approach would be necessary for these systems. 
We may also consider interacting topological defects in 
the patterned structures of driven (see, e.g., ]4l]| ) systems 
or of systems with competing interactions (see, e.g., [§)■ 
To apply our approach, the long-wavelength fluctuations 
(//), the interaction (n), the mobility T](l), and the nature 
of any noise-driven transport must be identified. 



VI. LEVY SUPER-DIFFUSION 

Systems with long-range Levy super-diffusion |Q have 
been used to model stirred reaction-diffusion systems. 
While there are no long-range interactions per se in these 
systems, super-diffusion enhances the reaction rate in a 
manner qualitatively similar to long-range interactions. 
Indeed, our methods can be applied to this case and agree 
with the results of Zumofen et al. ]l7j] for the late-time 
evolution from uncorrelated initial conditions. We also 
obtain additional information about the reaction zone 
structure and interface profiles for segregated morpholo- 
gies. For simplicity we only consider uncorrelated initial 
conditions, with p = d/2. 

In the discrete formulation of Levy flight, every particle 
hops a random distance r along a random lattice direc- 
tion with probability distribution P(r) oc r -1-7 — an- 
nihilating any anti-particles it encounters along the way. 
We impose 1 < 7 < 2 so that the hop distribution has a 
finite first moment and is normalizable. The equivalent 
continuum dynamics is dtPk — ~D\k\" 1 pk, where 7 > 1 is 
required. 

In time t, taken to be large, a particle randomly Levy 
walks a distance R ~ i 1 / 7 . If 7 > d, the volume bound- 
ing this walk R d <C t and the walk is recurrent, or space- 
filling. Conversely, if 7 < d the walk is sparse, and the 
particle explores a volume proportional to t using its fi- 
nite capture radius r c . Paralleling Eqn. (0), the anni- 
hilation time in a well-mixed region of the system scales 
as 



tl 



AB 

H 

AB 



J<d, 
7 > d, 



(14) 



where the typical particle-antiparticle separation in the 
region is Iab- From the Levy flight, the net distance that 

1 /'Y 

the particles move before annihilating is W ~ r L . For 
space-filling walks (7 > d) W ~ Iab, while for sparse 

walks (7 < d) W ~ l A ll > l AB - 

For a segregated morphology, there is a typical flux- 
density leaving domains. The flux-density is most eas- 
ily obtained in the discrete formulation, where it is sim- 
ply the number of particles that pass through a site in 
each time-step. From the Levy flight distribution, a 



fraction p ~ l/x 7 of particles contributes to the cur- 
rent from a given distance x away. Only particles along 
lattice directions contribute in a single time-step and, 
since 7 > 1, the flux is dominated by particles nearby 
compared to the domain scale. Hence, our current is 
J(r) ~ J dy[p(r + y) — p(r — y)]/y 7 a distance r from 
a domain edge, where W <C r <C L. Using the domain 
profile from Eqn. (|J), and requiring that J(r) approaches 
a constant near the domain edge we obtain a non-linear 
profile exponent, a = 7 — 1, and Jl ~ pL 1 1 for the 
particle current. 

For segregated systems, we equate the flux out of do- 
mains d t p ~ JL d ~ x /L d ~ Iaa +2 ^ to tne t™ e derivative 
of the charge density 1~aaI^ *° obtain Iaa ~ t 1 ^ 2 ^ and 
p ~ £-< i /( 2 i') Thg domain size L ~ i 1 / 7 . We can also 
obtain the rate of particle annihilation from the reaction 
rate within the reaction zones, dtp ~ WL d ~ x /[L d l AB TL], 
and compare to obtain 



AB 



£(d+2 7 -2)/[2d(2 7 -l)] 7 < d, 

t (d+2 7 -2)/[2 7 ( 7 +d-l)] 1> d. 



(15) 



Since particles are randomly Levy-walking before anni- 
hilation, the reaction zone width is W ~ Iab for d < 7, 
and W ~ l d /g for d > 7. 

The system has a segregated morphology, with W < 
L ~ l AA , for 7 > d/2. For 7 < d/2 the system is 
mixed, and the time-scale is set by the annihilation rate 
t l - t, so that Iaa ~ and p~l/t. At 7 = d/2, the 
marginal case between segregated and mixed morpholo- 
gies, we have W ~ L and Iab ~ Iaa- Our results agree 



with those of Zumofen et al. |17|. We also obtain a, Iab, 
and W characterizing the domain profile and reaction 
zone in segregated systems. 

It is interesting that while our results are qualita- 
tively similar to those obtained with long-range interac- 
tions, there is no effective interaction n that recovers the 
growth-exponents for a given Levy exponent 7. This is 
in contrast to the nonequilibrium steady-state properties 
of the kinetic Ising model, where Levy flights generate 
effective long-range interactions Q|. 

We may be able to apply these results to the motion 
of an single charged particle in a quenched random po- 
tential. This motion is sub-diffusive, with 7 > 2 and 
R ~ i 1 / 7 , for potentials with sufficiently long-range cor- 
relations f45|] . One loop RG calculations by Deem and 
Park [H for potentials caused by quenched Coulombic 
ions in d = 2 (n = 1), matched to the homogeneous evo- 
lution without any long-range interactions, indicate that 
sub-diffusive behavior alone may be sufficient to describe 
the evolution of the quenched system. If so, then our 
Levy-flight results may be directly applied with the ap- 
propriate 7. However significant questions remain, such 
as the appropriate microscopic annihilation mechanisms 
for oppositely charged particles in a quenched random 
potential. 







VII. LONG-RANGE CUTOFFS 

It is interesting to explore quenched long-range systems 
numerically, but the computational burden can be large. 
Long-range cutoffs L cut (t), usually on the order of the 
interparticle spacing, accelerate a simulation, but often 
at the cost of changing the physics. 

Long-range cutoffs contribute in two places. The first 
is through the force integral driving particles fluxes, F(r) 
in Eqn. (|l0|). Any L cut -C L(t) replaces the upper-cutoff 
of the integral and change F(r) for the cases where the 
integral is not dominated by short scales, i.e. for n < d + 
1. If F(r) is changed, the domain profile exponents a are 
also affected. The second place cutoffs enter is through 
the annihilation time within the reaction zone, Eqn. (12). 
If the cutoff is smaller than the reaction-zone spacing, 
L C ut <C Iab (t) i then the annihilation dynamics of Tp is 
changed. Ballistic annihilation times tb is qualitatively 
changed whenever L cut <C L(t). Amplitudes can also be 
affected by cutoffs, even when growth exponents are not. 

As an illustration, we consider limiting interactions to 
nearest- neighbors. Consider uncorrelated initial condi- 
tions with no diffusion, Fig. || a). For nearest-neighbor 
interactions, we use L cut ~ lAA( r ) ~ p( r )~ i the local 
particle spacing. Putting L cut instead of L in Eqn. fll(i| ) 
leads to the same results as n > d+1 in Table | for all n. 
In particular, the reaction-zone profile is characterized by 
a = d/(n+d— 1). The nearest- neighbor cutoff is of order 
Iab within the reaction zone so that local force-driven an- 
nihilation dynamics, with tf, is qualitatively unchanged. 
The system size cannot enter the force integral, so there 
is no restriction on the interaction to n > d/2. For a 
segregated morphology regime [i] applies, while in the 
mixed regime [Hi] does, though with W ~ Iab and anni- 
hilation through Tp. Comparing the particle annihilation 
rate dtp ~ J/L to the maximum rate supported by reac- 
tion zones dtp m ax ~ p/t f (Iaa) ~ 1/Iaa +1 , we nnd tnat 
segregation occurs for all d < 3. In summary, nearest 
neighbor interactions only leave the evolution qualita- 
tively unchanged for n > d + 1 . A numerical test of the 
effects of nearest-neighbor interactions in d = 1 is made 
in the next section. 



VIII. NUMERICAL INVESTIGATION 

We present some numerical results on the non- 
trivial domain profile exponent a in uncorrelated one- 
dimensional systems without diffusion. We also consider 
the effect of a cutoff L cut limiting interactions to nearest- 
neighbors. More results will be presented in future work 
& 

In d = 1, we expect a = (n — l)/2 for 1 < n < 2 
(regime [ii ]) and a = 1/n for n > 2 (regime [i]). We 
studied systems with 8000 particles with n = 3/2 [375 
samples] and n — 2 [300 samples]. We expect a = 1/4 
and 1/2 respectively. We avoided finite-size effects by 



comparisons with size 4000 systems. We also considered 
systems with only nearest-neighbor interactions, which 
are e xpec ted to modify the domain profile for n < 2, see 



VII 



Sec. 
for all n 



In particular, 
< 2. For n 



we expect regime 
3/2, we expect a 



[i\ to apply 
= 2/3. We 

studied systems with 128000 particles with n = 3/2 [200 
samples] and n = 2 [53 samples], and checked finite-size 
effects with systems of 32000 particles. We found results 
consistent with expectations. We show our results for the 



domain profiles in Fig. ^|. 



IX. CONCLUSIONS 

By considering the annihilation dynamics in well mixed 
regions of a charged particle system, and balancing the 
annihilation against currents driven by charge inhomo- 
geneities left over from the initial conditions, we self- 
consistcntly determine the morphology and evolution of 
quenched charged-particle systems with long-range inter- 
actions. We also contribute a visceral description of the 
dynamics. We characterize the system with the scale of 
domains L and reaction-zones W, and the particle spac- 
ings within domains Iaa and reaction-zones Iab- For 
mixed systems L ~ Iaa ~ Iab- Our results are sum- 
marized in Fig. and in the three tables. Our primary 
assumption is that the lengths we have used are suffi- 
cient to characterize the evolving system. The scaling 
form for the domain profile, Eqn. (||), follows from this 
assumption. 

The results of this paper will hopefully inspire more 
formal derivations and numerical tests, as well as exper- 
imental tests in electronic systems. Comparisons with 
existing treatments is encouraging, particularly agree- 
ment with the field-theory approach of Lee and Cardy 
Q for reaction-diffusion systems, with the hydrodynamic 
treatment of Ginzburg, Radzihovsky, and Clark |L2| for 
uncorrelated initial conditions with diffusion and long- 
range interactions, and with the Levy super-diffusion re- 
sults of Zumofen et al. p7[ . There are also many new 
results pertaining to the reaction-zones, domain profiles, 
and systems with equilibrated high-temperature initial 
conditions. 

Electronic systems (n = 2) in two-dimensions should 
provide an experimental test for our results. Asymptot- 
ically, we predict a segregated morphology with domain 
size L ~ t 1 / 2 and average density p ~ £~ 3 / 4 for photo- 
excited quantum-well or QHE systems (high-temperature 
equilibrium initial conditions, regime [viii' / viii"]) . This 
will be explored at more length in a separate publica- 
tion ]35| ]. If initial conditions with uncorrelated initial 
conditions can be manufactured, then we expect regime 
[iv' /iv" /ii /ii 5 ) to apply: a segregated morphology with 
L ~ i 1 / 2 and p ~ t~ 1//2 . These contrast dramatically 
with the mixed morphology and density decay p ~ t^ 1 
in three-dimensional electronic systems. 
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FIG. 1. Charge density fluctuations 8p ~ L~ M for equi- 
librated high-temperature initial conditions. In the shaded 
region the initial conditions are uncorrelated at large scales; 
in the unshaded region interactions reduce large-scale charge 
fluctuations. 




FIG. 2. The asymptotically dominant flux starting from 
uncorrelated initial conditions (/i = d/2). When Jf domi- 
nates, see Tab. |. when Jd dominates see Eqn. For equi- 
librated high-temperature initial conditions Jd always domi- 
nates at late times. 



FIG. 3. Annihilation time t(1) vs size / of a parti- 
cle/anti-particle pair placed in periodic boxes of size I in 
d = 3, where the particles are initially 1/4 apart. They evolve 
according to Eqn. (^) with a central force / = l/r n , mobility 
T) = 1, diffusion constant D = 1/4, and time-step St = 1. The 
dashed curves correspond to, from bottom to top, n — 1/2, 1, 
3/2, 2, 3, and 4. Triangles correspond to purely diffusive mo- 
tion. The short solid lines indicate the expected asymptotic 
behavior, with t(1) ~ l 3 ^ 2 , I 2 , and I 3 (twice), from bottom to 
top respectively. Diffusion, with td ~ I 3 , dominates asymp- 
totically for n > 1. 




FIG. 4. Dominant annihilation times for force and noise, 
with uncorrelated initial conditions. Diffusive annihilation 
(td) dominates in the shaded region, while ballistic annihi- 
lation (jb ~ I'abF' 1 ^ 2 ) dominates in the clear region. The 
regions [ii 1 ] and [ii 2 ] correspond to Fig. ^. For equilibrium 
high-temperature initial conditions, when noise is present, td 
always dominates. For force-only annihilation, irrespective of 
initial conditions, tb ~ l d AB P i ' d ~ n ~ 1 ^ n dominates for n < d 
while tf ~ ZTt dominates for n > d. 
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FIG. 5. Different regimes of density evolution. Blank regions indicate a segregated morphology, shaded a mixed morphology, 
while black indicates excluded regimes that have no thermod yna mic limit. Dashed lines separate regimes, numbered with Roman 
numerals, with different reaction-zone structure (see Table |lll|). The figures correspond to: a) uncorrelated initial conditions, 
force only, b) uncorrelated initial conditions, diffusion and force, c) enlargement of b), with the [ii 1 ] to [ii 5 ] subregions numbered, 
d) equilibrated high-temperature initial conditions, force only, and e) equilibrated initial conditions, noise and force. 
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profiles: a = 1/4, 2/3, and 1/2. The crossover at small x is due to the reaction zone, and moves to smaller scaled distance at 
later times. 
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r (d- n )/2^(d- n -2ti)/2 


j^d — n — 2jj, 


n > d+ 1 


d/(n + d- 1) 


(L/r) d/(n+d_1) L _(n+1)/2 


£-(n + d+l)/2 



TABLE I. Domain profile exponents, a, coarse-grained field a distance r<i from the domain edge, F(r), and force-driven 
flux, Jf, for different interaction exponents n. When diffusive fluxes dominate a = 1. 
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TABLE II. Growth laws of Iaa, L, Iab, and J for the different labeled regimes in Fig. g. 
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TABLE III. Growth laws of W, t, and h for the different labeled regimes in Fig. [5|. Reaction-zone widths I4 7 with the same 
scaling as Iab are indicated. The dominant mechanisms determining W and r are also indicated, where "D" is diffusive, "F" 
is force driven, and "B" is ballistic. Note that when the morphology is mixed, r ~ t and W indicates the scale of remaining 
charge fluctuations. 
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